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SUMMARY 


INTRODUCTION 

When  a  force  is  impressed  upon  a  body  immersed  in  a  fluid  continuum, 
the  resultant  motion  of  the  body  is  determined  not  by  the  impressed  force 
alone,  as  it  would  be  in  vacuo,  but  by  the  impressed  force  and  the  forces 
arising  from  the  body -"Fluid  interaction. 

If  the  body  moves  linearly  with  constant  velocity  through  an  other¬ 
wise  undisturbed  homogeneous,  incompressible  fluid  with  sufficiently  re¬ 
mote  boundaries,  or  if  there  is  a  uniform  (constant  velocity)  flow  of  the 
fluid  past  the  body,  the  interaction  force  can  be  described  as  a  resist¬ 
ance,  or  drag.  However,  the  body-fluid  interaction  for  unsteady  motions 
is  not  adequately  described  by  the  drag  for  uniform  linear  motion.  Cor¬ 
responding  relationships  must  be  developed  for  the  accelerating  body.  An 
attempt  to  describe  the  acceleration-induced  resistance  leads  to  the  con¬ 
cept  of  virtual  mass. 

Motion  of  a  body  in  a  fluid  medium  results  in  the  interdependent 
phenomena  of  boundary  layer,  wake,  and  mass  transport.  For  steady  motion, 
a  total,  or  composite,  drag,  which  includes  deformation  drag,  friction 
drag,  and  pressure  drag,  is  defined,  expressing  the  resistance  to  motion 
as  a  quadratic  law  containing  an  empirical  nondimensional  drag  coeffi¬ 
cient.  The  total  drag  depends  on  the  velocity,  body  geometry,  and  phys¬ 
ical  properties  of  the  fluid,  the  parameters  controlling  the  wake  and 
boundary  layer.  Although  a  net  fluid  mass  transport  occurs,  it  does  not 
contribute  to  the  drag,*  Only  when  acceleration  takes  place  does  the 
effect  of  transported  fluid  mass  become  a  factor  in  resistance. 

When  a  body  is  accelerated  in  a  fluid,  an  inertial  force  occurs 
because  a  mass  of  fluid  is  also  accelerated.  In  otner  words,  the  force 
which  acts  to  impart  kinetic  energy  to  the  body,  through  the  body-fluid 
interaction,  imparts  kinetic  energy  to  the  surrounding  fluid,  as  well. 

The  total  mass  being  accelerated,  or  apparently  being  accelerated, 
is  called  the  virtual  mass.  Tnerefore,  the  virtual  mass  is  the  sum  of 
the  physical  mass  of  the  body  and  the  hydrodynamic  mass  (often  termed 
"added"  mass  or  "induced"  mass) .  The  hydrodynamic  mass  corresponds  to  a 
volume  of  fluid  defined  conventionally  in  simple  geometrical  terms.  The 
simplified  concept  of  an  entrained  fluid  mass  moving  with  the  body  rep¬ 
resents  a  complex  phenomenon  of  body-fluid  interaction  which  plays  an 
important  part  in  hydrodynamics. 

Sonobuoy  systems  consist  of  essentially  three  elements:  a  surface 
float,  a  hydrophone,  and  an  intermediate  suspension.  As  the  surface 
float  is  required  to  maintain  an  exposed  antenna,  it  is  designed  to 
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follow  the  ocean  surface  with  a  constant  freeboard.  The  hydrophone,  on 
the  other  hand,  is  susceptible  to  adverse  dynamic  effects  and  is  required 
to  be  virtually  motionless.2.3  The  function  of  the  intermediate  suspen¬ 
sion  is  to  prevent  the  transmission  of  motion  of  one  to  the  other. 

A  hydrophone  suspension  invariably  takes  the  form  of  a  simple  spring- 
mass  system  or  series  of  spring -mass  systems.  It  is  so  designed  that  its 
natural  frequency  is  below  the  range  of  ocean  wave  frequencies,  and  there¬ 
fore,  the  motion  amplitude  of  the  hydrophone  is  attenuated  from  that  of 
the  surface  float.  The  smaller  the  natural  frequency  of  the  suspension, 
the  more  effective  the  suspension  is. 

Three  parameters  affect  the  natural  frequency  -  mass,  spring  con¬ 
stant,  and  damping  coefficient.  For  a  low  natural  frequency,  mass  arid 
damping  coefficient  should  be  large,  and  spring  constant  small.  The  jug¬ 
gling  of  these  parameters  for  optimum  efficiency  of  motion  isolation  is 
influenced  to  a  certain  extent  by  size,  weight,  and  packaging  considera¬ 
tions  . 

Large  physical  mass  means  a  heavy  sonobuoy,  but  beyond  this,  a 
weighty  lower  unit  limits  the  use  of  a  nighly  elastic  compliance  in  the 
suspension.  To  overcome  this  limitation  on  the  compliance,  buoyancy  can 
be  utilized  to  reduce  weight  while  retaining  mass.  The  use  of  buoyancy, 
however,  vastly  increases  size,  possibly,  beyond  packageability.  One 
means  of  acquiring  mass  without  significant  weight  or  size  problems  is 
the  utilization  of  hydrodynamic  mass. 

Practical  utilization  of  hydrodynamic  mass  for  the  design  of  a  sono- 
buoy  suspension  cannot  be  implemented  until  the  phenomenon  itself  has 
been  sufficiently  understood  and  described  for  design  criteria  to  be  de¬ 
termined.  These  criteria  have  not  heretofore  been  available.  The  task 
of  ascertaining  these  criteria  was  undertaken  at  NAVAIRDEVCEN  under  In¬ 
dependent  Research  Task  No.  ZR011-01-01,  Work  Unit  No.  ES-1-01, 

RESULTS  AND  CONCLUSIONS 

1.  The  hydrodynamic  mass  and  acceleration  resistance  effects  of  oscil¬ 
lating  disks,  cylinders,  and  tandem  disks  have  been  determined  through 
measurement  of  the  frequency  and  amplitude  of  motion  at  resonance  of  a 
..spring-mass  system. 

2.  The  hydrodynamic  effects  on  oscillating  disks  (and  other  bodies) 
can  he  expressed  in  terms  of  coefficients  of  hydrodynamic  inertia,  drag, 
and  total  resistance,  and  can  be  correlated  by  the  acceleration  modulus 
v2/ad,  which  is  a  measure  of  motion  amplitude,  although  no  correlation 
with  Reynolds  number  is  apparent. 


2,3  See  References  .■  nage  25. 
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3.  The  hydrodynamic  coefficients  describing  a  disk  in  harmonic  accel¬ 
eration  can  be  written  as: 


Ci  =  m^/pd3  *  1.2  (v2/ad)1/2 


Cq  =  FD/(pAv2/2)  =2.2  (vVad)-1/2 


CT  =  FT/(pAv2/2)  =  5.2  (v2/ad)-1/2 
within  limits  of  motion  approximated  by  0.08  v2/ad  3. 

4.  The  hydrodynamic  inertial  effect  is  associated  with  the  vortex  for¬ 
mation  in  the  wake  of  the  accelerated  body,  and  an  anomaly  occurs  wh^n 
the  vortex  is  shed,  causing  unstable  motion  of  the  oscillating  body. 

5.  Variation  of  hydrodynamic  inertial  effects  of  oscillating  cylinders 
with  length-to-diameter  ratio  and  of  oscillating  tandem  disks  with  sepa- 
tation-to-diameter  ratio  has  been  determined. 

6.  The  data  resulting  from  this  investigation  can  be  applied  directly 
to  sonobuoy  suspension  systems,  facilitating  the  analysis  and  design  of 
these  systems. 
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NOMENCLATURE 

A  =  Projected  surface  area  normal  to  motion 
a  -  Acceleration 

B  =  Volume  characteristic  of  a  body 

Cp  =  Drag  coefficient 

Cj  =  Hydrodynamic  inertia  coefficient 

Cx  -  Total  resistance  coefficient 

c  =  Damping  constant 

cc  =  Critical  damping  constant 

d  =  Diameter  of  body 

f  =  Frequency 

Fp  =  Drag  force 

Fj  =  Hydrodynamic  inertia  force 
Fj  =  Total  resistance  force 
g  =  Gravitational  acceleration 
k  =  Spring  constant 
i  =  Cylinder  length 
M  =  Magnification,  y/y0 

m  =  Virtual  mass,  m^  +  mp 

=  Hydrodynamic  mass 
mp  =  Physical  mass 
Np  =  Froude  number 
Npe  =  Reynolds  number 
St  =  Strouhal  number 

-  viii  - 
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NOMENCLATURE  (continued) 

s  =  Separation  distance  between  tandem  disks 

T0  *  Oscillation  period 

Ts  *  Vortex  shedding  period 

t  *  Time 

v  =  Velocity 

x  =  Linear  distance 

y  *  Amplitude  of  oscillatory  motion 

y0  =  Input  amplitude  to  spring-mass  system 

v  =  Kinematic  viscosity 

p  =  Fluid  density 

a)  =  Angular  frequency,  2irf 

*  Undamped  natural  frequency  of  spring-mass  system  in  air 
Unj  =  Damped  natural  frequency  of  spring -mass  system  in  fluid 
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BACKGROUND 

The  effect  of  acceleration  on  a  body  in  a  fluid  was  first  noted  in 
the  1776  work  of  DuBuat  with  spherical  pendulum  bobs  in  water;  then,  ir. 
the  early  nineteenth  century,  Bessel,  Green,  and  Stokes  laid  the  theoret¬ 
ical  foundation  for  the  concept  of  virtual  mass.  Extensive  bibliographies 
on  early  work  in  this  subject  have  been  presented  by  Torobir.  and  Gauvin*4 
and  Stelson.5 

The  theoretical  basis  for  hydrodynamic  mass  arises  out  of  potential 
flow  considerations  of  the  kinetic  energy  imparted  by  a  moving  body  to  an 
ideal,  incompressible  fluid  of  infinite  extent.  Milne-Thomson6  considered 
the  disk  as  a  special  case  of  the  planetary  ellipsoid.  From  the  complex 
potential,  the  kinetic  energy  was  determined  to  be  pd3v2/6,  for  a  circular 
disk  of  diameter,  d,  moving  perpendicularly  to  its  plane  with  velocity, 
v,  in  a  fluid  of  density,  p.  The  hydrodynamic  mass,  therefore,  is 


=  pd3/3  (i) 

Kanwal7  considered  oscillating  axially-symmetric  bodies  at  small 
Reynolds  numbers  to  determine  virtual  mass  arising  in  Stokes  flow.  The 
formula  derived  for  a  circular  disk  moving  broadside  to  the  stream  for 
small  Reynolds  numbers  exhibited  a  dependency  on  the  square  root  of  os¬ 
cillation  period.  A  similar  result  was  obtained  for  the  sphere,  agreeing 
with  the  parameter  presented  by  Carstens®  to  describe  spherical  particle 
motion  in  a  fluid. 

The  treatment  of  an  oscillating  disk  by  the  integration  of  pressures 
method  based  on  a  vibrating  circular  piston  in  an  infinite  baffle9  results 
in  a  series  solution  for  hydrodynamic  mass.  At  high  frequencies,  the 
value  of  hydrodynamic  mass  is  approximated  to  exhibit  a  dependence  on  the 
square  of  oscillation  period,  and  at  lower  frequencies,  mjj  approaches 
pd3/3,  the  same  value  obtained  in  equation  (1).  The  assumption  of  this 
method  is  that  the  amplitude  of  oscillation  be  extremely  small,  since  no 
account  is  taken  of  boundaiy  layer  formation. 

The  conditions  of  an  incompressible,  irrotational  ideal  fluid  under 
which  the  results  of  potential  flow  theory  apply  are  difficult  to  achieve 
with  real  fluids  in  any  but  very  limited  situations.  For  a  circular  disk 
in  flowing  water,  Stanton  and  Marshall10  found  a  vortex-ring  to  have 
formed  at  a  Reynolds  number  of  5,  and  vortex  shedding  to  have  begun  at 
Nge  =  195.  Simmons  and  Dewey,11  Willmarth,  et  al,  and  others  have  re¬ 
ported  the  critical  Reynolds  number  for  vortex-related  instability-onset 
in  disk  motion  as  100.  Once  vortex  formation  occurs,  the  theoretical 
solutions  have  no  great  practical  importance. 


5-12  See  Refer ances,  page  2b. 
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Small  amplitude  oscillations  do  approximate  the  conditions  of  the 
theory,  and  the  way  in  which  the  hydrodynamic  mass  behaves  in  this  area 
is  of  interest.  Stelson5  reported  on  experiments  with  bodies  oscillating 
in  fluid  with  amplitudes  on  the  order  of  0.01  inch  at  frequencies  between 
14  and  35  Hz.  For  2-inch  diameter  disks,  the  values  of  hydrodynamic  mass 
obtained  were  in  good  agreement  with  the  theoretical  value,  being  slightly 
larger  when  disk  thickness  was  increased.  Similarly  small  amplitudes 
(0.002  to  0.011  diameter)  for  disks  vibrating  in  water  at  10.31  Hz  were 
employed  by  Bramig13  who  reported  values  of  hydrodynamic  mass  below  the 
theoretical  value  but  approaching  it  as  the  amplitude  decreased .  Carstens6 
noted  the  effect  of  increasing  amplitude  of  oscillation  on  submerged 
spheres  to  be  a  reduction  of  hydrodynamic  mass  from  the  theoretical  value. 

The  indication,  in  the  works  of  both  Bramig  and  Carstens,  that  hy¬ 
drodynamic  mass  decreases  as  the  amplitude  of  oscillation  increases,  is 
misleading,  because  larger  amplitude  oscillation  experiments  performed 
more  recently  denote  the  opposite  trend.  This  apparent  contradiction 
probably  occurs  because  of  the  different  flow  regimes  in  which  the  exper¬ 
iments  took  place. 

The  conditions  of  turbulent  flow,  which  accompany  any  appreciable 
notion  of  a  bluff  body  in  a  real  fluid,  have  been  encountered  in  experi¬ 
ments  involving  both  linear  motion  and  oscillation.  Significant  results 
have  been  obtained  for  linear  motion  by  Iversen  and  Balent14  and  Luneau15 
for  disks,  Keim16  for  cylinders,  and  Bugliarello17  and  Lunnon18  for 
spheres.  Oscillatory  motion  has  been  investigated  for  cylinders  and  flat 
plates  by  Keulegan  and  Carpenter18  and  McNown  and  Keulegan,20  for  cylin-  jj 

ders  and  spheres  by  Crooke,21  and  for  various  bodies  by  Patton.22  Many  ? 

of  these  investigations  have  been  summarized  by  Torobin  and  Gauvin4  and 
by  Wiegel.23  The  hydrodynamic  mass  determined  in  these  investigations 
for  several  bluff  bodies  in  two  distinct  modes  of  accelerated  motion  was 
larger  than  the  theoretical  value  for  potential  flow  and  showed  a  distinct 
anomally  when  vortex  shedding  was  initiated.  Further  reference  to  these 
investigations  is  made  in  the  subsequent  discussion  of  the  results  of  the 
present  study. 


13-23  See  References ,  pages  25  and  26. 
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DISCUSSION 
DETERMINATION  OF  HYDRODYNAMIC  MASS 

The  experimental  technique  used  to  measure  the  hydrodynamic  mass  of 
bodies  in  high  amplitude  oscillation  was  the  natural  frequency  method. 
This  consisted  of  measuring  the  shift  in  resonant  frequency  and  the  am¬ 
plification  of  body  excursion  for  the  spring-supported  body  in  water  and 
air.  In  this  manner,  the  effect  of  fluid  damping  was  determined  at  the 
same  time  as  hydrodynamic  mass. 

Apparatus 

The  experiments  were  undertaken 
using  1/16-inch  thick  aluminum  disks 
of  2,  4,  6,  and  12-inch  diameters, 
and  a  12-inch  diameter  fiberglas 
disks  of  the  same  thickness.  The 
disk  was  mounted  orthogonally  on 
the  end  of  a  1/4-inch  diameter  rigid 
rod,  attached  at  its  center,  and  the 
rod  was  suspended  from  a  helical 
spring.  A  variable  speed  drive, 
with  adjustable  input  (forcing)  am¬ 
plitude  from  1.5  to  5.5  inches,  was 
used  to  oscillate  the  top  end  of  the 
spring  through  a  range  of  frequen¬ 
cies.  The  experimental  setup  is  il¬ 
lustrated  in  figure  1. 

The  springs  were  fabricated  of 
piano  wire  wound  in  close  coils. 

The  helical  springs  were  designed  to 
propriate  masses,  the  springs  exhibited  internal  damping  ratios  c/cc 
<  0.005.  The  springs  were  determined  to  be  linear  and  the  spring  con¬ 
stants  were  invariable. 

.  The  water  tank  was  cylindrical,  having  a  depth  of  16  feet  and  a  di¬ 
ameter  of  20  feet.  The  disks  were  submerged  in  tap  water  to  a  mean  depth 
of  4  feet.  Experiments  were  spheres  by  Stelson5  indicated  2  diameter 
submergence  produced  negligible  free  surface  effects.  Murtha24  places 
this  depth  limit  at  2.5  diameters  for  cylinders;  Patton,22  at  5  diameters 
for  spheres;  *nd  Waugh  and  Ellis,25  at  2.5  diameters  for  spheres.  Of 
these  investigators,  only  Patton  experienced  high  amplitude  oscillation. 
Sufficient  submergence  depth  was  considered  tc  have  been  obtained  for  the 
present  experiments  so  long  as  the  surface  was  not  disturbed  by  the  disk 
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22-25  See  References,  page  26. 


-  3  - 


motion.  When  the  disk  was  too  close  to  the  water  surface,  a  definite 
rippling  accompanied  each  motion  cycle,  and  the  resonant  peak  shifted  to 
a  higher  frequency  than  that  observed  when  the  disk  was  at  sufficient 
depth.  Good  agreement  was  obtained  between  the  data  from  disks  in  this 
large  tank  and  those  from  similar  preliminary  tests  in  a  2-  by  2-foot 
square  tank,  having  a  depth  of  3.5  feet,  with  disks  up  to  6-inch  diameter. 

Additional  tests  wore  conducted  using  cylinders  of  2.5-inch  diameter 
and  length-to-diameter  ratios  of  one  and  two.  These  were  oscillated  ax¬ 
ially.  Two  similar  disks  in  tandem  were  also  employed. 

Procedure 

The  spring-mass  system  was  assumed  to  be  adequately  described  by  the 
linear  single-degree-of-freedom  equation: 

my  +  cy  +  ky  =  F0  cos  tot,  (2) 

where  m  the  virtual  mass  (the  sum  of  physical  mass  nip  and  hydrodynamic 
mass  mh)  and  c  the  damping  coefficient  for  the  disk  in  water  are  dependent 
upon  the  system  parameters.  The  peak  amplitude  and  the  frequency  at  which 
it  occurred  were  determined  experimentally.  This  resonant  frequency  was 
considered  the  damped  natural  frequency  for  the  system  in  watt  *  and 
comparison  with  the  undamped  natural  frequency  in  air  yield  equations 
for  the  hydrodynamic  mass: 

mh  -  mp  ([(2  *  Q)/6]  K/^dl2  -  U,  (3} 

and  coefficient  of  damping: 

c  =  (mpto'/tond)  /(4  -  Q)(2  ♦  Q)/9,  (4) 


where  Q  (equation  (A- 14) ;  appendix  A)  depends  solely  upon  the  magnifica¬ 
tion  factor  of  the  motion  amplitude.  Derivation  of  equations  (3)  and  (4) 
are  shown  in  appendix  A. 

HYDRODYNAMIC  MASS  AND  DAMPING  OF  OSCILLATING  DISKS 

The  drag  of  a  bluff  body  in  uniform  flow  is: 

Fd  =  CDpAv2/2,  (5) 

where  Cq  is  the  nondimcnsional  drag  coefficien.,  p  is  the  fluid  density, 
a  is  the  projected  surface  area  of  the  body  normal  to  the  direction  of 
flow,  and  v  is  the  velocity  of  the  body  relative  to  the  fluid.  For  a 
body  of  a  given  shape,  the  drag  is  described  by  equation  (5)  and  a  plot 
of  Cp  versus  Reynolds  number,  as  determined  by  experiment. 
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Resistance  to  acceleration  may  be  considered  as  a  hydrodynamic  iner¬ 
tial  force: 

Fj  =  raf,a  =  CjpBa,  (6) 

where  B  is  a  volume  characteristic  of  the  body,  a  is  the  body  accelera¬ 
tion,  and  Cj,  like  Cq,  may  be  expected  to  be  dependent  upon  the  condition 
of  body  motion. 

Because  most  investigations  of  hydrodynamic  acceleration  effects 
consider  linear  motion  where  the  inertial  and  drag  forces  act  in  concert, 
a  variation  of  the  resistance  force  equations  has  arisen: 

F"p  -  CqpAv2/2  ♦  CjpBa  =  CjpAv2/2,  (7) 

where  Fj  is  a  total  resistance  force,  and  C-p  is  a  total  resistance  coef¬ 
ficient.  In  oscillatory  motion  where  velocity  and  acceleration  have  a 
tt/2  phase  difference,  the  concept  of  a  simple  summation  of  drag  and  hydro- 
dynamic  inertia  as  if  they  were  a  single  drag  force  would  seem  of  dubious 
merit,  considering  the  time-sequence  of  harmonic  motion.  However,  the 
values  of  hydrodynamic  inertia  and  drag  obtained  in  this  investigation 
resulted  from  evaluation  of  the  motion  cycle  of  the  oscillating  body  and 
may  be  said  to  represent  steady-state  equivalents  rather  than  instanta¬ 
neous  values.  As  both  the  drag  and  inertial  forces  are  acting  in  opposi¬ 
tion  to  the  body  motion  throughout  the  oscillation  cycle,  they  can,  in  a 
sense,  be  summed  into  the  single  resistance  force  in  equation  (7).  Values 
for  Cj  have  been  calculated  for  the  oscillating  disk  in  this  investigation 
for  comparison  with  linear  motion  data  reported  in  the  literature;  how¬ 
ever,  this  procedure  must  be  accompanied  by  appropriate  notation  of  the 
phase  discrepancy  involved  in  placing  the  oscillatory  data  on  the  same 
basis  as  the  linear  data. 

Because  the  linear  single-degree-of-freedom  equation  was  used  for 
the  calculation  of  hydrodynamic  mass  and  damping  coefficient  c,  a  drag 
coefficient  corresponding  to  the  turbulent  damping  was  determined  by 
equating  the  damping  force  cv  to  a  drag  force: 


CV  =  CqpAv2/2, 

(8) 

defining  the  drag  coefficient  as: 

Cq  *  2c/pAv. 

(9) 

This  value  was  employed  in  equation  (7)  for  the  calculation  of  Cq*. 

According  to  the  theory  of  potential  flow,  the  hydrodynamic  mass  is 
the  mass  of  a  volume  of  fluid  equivalent  to  63.7  percent  of  the  volume  of 
a  sphere  with  the  same  diameter  as  the  disk.  This  can  be  expressed  more 
simply  as  pd3/3.  Therefore,  the  volume  B  in  equation  (6)  is  chosen  as  d3 
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rather  than  the  volume  of  a  sphere,  as  the  spherical  volume  has  no  partic¬ 
ular  physical  significance  for  the  disk.  The  inertia  coefficient  Cj  for 
potential  flow  is  simply  one-third,  and  in  general: 

Ci  =  mh/pd3  (10) 

Correlating  Modulus 

Dimensional  analysis  yields  a  number  of  nondimens ional  groups  de¬ 
scriptive  of  the  body-fluid  interaction.  Those  groups  that  are  possibly 
meaningful  to  this  investigation  include: 

1.  Reynolds  number,  NRe  *  vd/v; 

2.  Froude  number,  Np  *  v//g cT; 

3.  Acceleration  modulus,  ad/v2; 

4.  Acceleration  rate  of  change  modulus,  ad2/v3;  and 

5.  Time  mcdulus,  vt/d. 

One  or  several  of  these  dimensionless  groups  may  be  useful  in  correlating 
the  data  in  the  form  of  the  nondimensions  1  coefficients  Cj,  CD,  and  CT. 

Free  surface  effects  were  found  to  be  negligible  in  this  investiga¬ 
tion,  as  the  test  body  was  submerged  to  sufficient  depth.  Therefore,  the 
Froude  number  was  not  an  influence  upon  the  resistance  to  acceleration. 

Reynolds  number  reflects  the  influence  of  viscous  forces  and  is  im¬ 
portant  in  the  correlation  of  drag  coefficients  for  uniform  flow.  The 
drag  coefficients  for  disks  and  thin  plates  in  constant  velocity  motion 
are  shown  in  figure  2.  Using  the  maximum  harmonic  velocity  for  calcula¬ 
tion  of  the  Reynolds  number,  the  values  of  Cj,  Cq,  and  Op  for  the  oscil¬ 
lating  disk  are  plotted  in  figures  3,  4,  and  5  against  Reynolds  number. 

No  distinct  correlation  occurs.  This  agrees  with  the  results  of  most 
investigators . 

Torobin  and  Gauvin,4  in  summarizing  tbs  literature  of  accelerated 
motion  effects,  noted  a  trend  for  the  effect  of  acceleration  on  drag 
coefficient  to  be  reduced  as  the  Reynolds  number,  and  therefore,  the 
turbulent  component  of  the  wake  structure,  was  increased.  There  is  not 
sufficient  evidence  in  the  data  of  this  'nvestigation  to  substantiate 
that  reported  trend. 

The  acceleration  modulus,  first  employed  by  Iversen  and  Balent,14 
has  proved  valuable  in  correlating  the  data  of  acceleration-induced  re¬ 
sistance.  The  modulus  ad/v2  has  a  significance  that  has  been  overlooked 
or  ignored  by  investigators  in  this  area.  To  discuss  this  point,  it  is 
convenient  to  define  the  reciprocal  acceleration  modulus  v2/ad. 
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FIGURE  2  -  Drag  Coefficient  Versus  Reynolds  Number  for  Disks  and  Thin 
Plates  in  Steady  Flow 


FIGURE  3  -  Hydrodynamic  Inertia  Coefficient  Versus  Reynolds  Number  for 
Circular  Disks  in  Harmonic  Acceleration 
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Considering  the  motion  of  a  body  accelerated  from  rest  at  a  constant 
acceleration,  then: 

a  *  *x  *  const 

v  «  i  »  at 

x  *  x  *  at2/2.  (11) 

The  Modulus  reduces  to  a  Measure  of  the  distance  travelled: 

v2/ad  *  2x/d.  (12) 

It  zb--,  be  shown  that  similar  results,  with  differences  in  proportionality 
constants,  can  be  obtained  for  linear  Motion  associated  with  practically 
any  continuously  acting  none tonic  acceleration. 

In  the  case  of  oscillatory  Motion,  naxiaun  values  are: 

a  *  y  =  ya2; 

v  *  y  *  y«; 


y  *  y;  'is) 

and  the  acceleration  Modulus  reduces  to: 

v2/ad  *  y/d.  (14) 

This  result  imparts  a  physical  significance  to  the  v2/ad  Modulus  and  in¬ 
dicates  the  relationship  between  linear  and  oscillatory  notions. 

The  rate  of  change  of  acceleration  modulus  for  oscillation  is  simply: 

v3/ad2  =  (y/d)2,  (15) 

and  can  be  considered  within  the  fraaework  of  the  modulus  v2/ad. 

The  ti«e  modulus  vt/d  has  been  used  by  Keulegan  and  Carpenter19  with 
success  in  correlating  the  data  of  plates  and  cylinders  in  oscillating 
flow.  For  oscillation,  it  is  apparent  that  time  modulus  and  acceleration 
modulus  are  closely  related: 

vt/d  *  v/fd  =  yw/fd  =  2vy/d.  (16) 

As  it  would  be  superfluous  to  employ  both  the  time  modulus  and  the 
acceleration  modulus  to  describe  the  motion,  the  succeeding  discussion 
considers  only  the  v2/ad  modulus. 
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Inertial,  Drag,  and  Resistance  Coefficients 

The  hydrodynamic  inertia  coefficient  Cj  for  oscillating  disks  is 
shown  in  figure  6  as  a  function  of  the  modulus  v2/ad.  Correlation  of  the 
sort  that  could  not  be  obtained  for  Reynolds  number  does  occur  for  the 
acceleration  modulus.  Using  the  method  of  least  squares,  an  approximate 
equation  for  the  inertia  coefficient  was  determined? 

Ci  =  1.2(v2/ad)  l/2  (17) 

for  0.08  <_  v2/ad  ^3.  For  values  of  v2/ad  <  0.08,  the  coefficient  Cj  = 
1/3,  the  theoretical  value.  When  v2/ad  >  3,  the  disk  motion  became  er¬ 
ratic  and  measurements  could  not  be  made. 

Iversen  and  Balent14  and  Luneau15  investigated  the  acceleration  ef¬ 
fects  on  circular  disks  in  linear  motion.  The  former  operated  with  vari¬ 
able  acceleration  (decreasing  with  time)  and  the  latter  with  constant 
acceleration.  Thei-  results  are  compared  with  those  of  the  present  in¬ 
vestigation  in  figuie  7.  It  should  be  noted  that  the  three  curves  rep¬ 
resent  different  relationships  between  v2/ad  and  y/d.  For  example, 
Luneau's  v2/ad  modulus  is  twice  the  distance  travelled  divided  by  disk 
diameter  [equation  (12)],  and  the  v2/ad  modulus  for  the  present  investi¬ 
gation  is  one-half  the  total  excursion  of  the  oscillating  disk  divided  by 
the  diameter  [equation  (14)].  If  these  two  curves  were  plotted  on  a 
scale  cf  equivalent  distance  travelled  in  diameters,  they  would  fall  much 
closer  together. 

When  the  damping  constant  c  of  equation  (2)  is  plotted  against  v2/ad, 
no  correlation  occurs  (figure  8).  However,  when  the  damping  constant  is 
converted  into  a  drag  coefficient  by  equation  (9),  correlation  with  v2/ad 
results,  as  shown  in  figure  9. 

Employing  equation  (7),  a  total  resistance  coefficient  C7  was  plotted 
in  figure  10.  A  least  squares  computation  lead  to  the  approximation: 

CT  =  5.2(v2/ad)-1/2,  (18) 

From  equations  (7),  (17),  and  (18): 

C3  =  2.2(v2/ad)'l/2.  (19) 

A  comparison  of  the  data  of  the  present  investigation  for  oscillating 
disks  with  the  lineal  motion  data  for  disks  from  Iversen  and  Balent  and 
Luneau  is  shown  in  figure  11.  The  conditions  of  potential  theory  for 
hydrodynamic  mass  and  of  steady  flow  drag  are  shown  as  limiting  valiuo 
after  the  technique  of  Burgliarello, 17  whose  data  on  spheres  followed 
similar  trends. 
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FIGURE  6  -  Hydrodynamic  Inertia  Coefficient  Versus  Acceleration  Modulus 
for  Oscillating  Disks 


FIGURE  7  -  Comparison  of  Data  for  Oscillating  Disks  and  Disks  in  Linear 
Acceleration 
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Disks  of  square,  hexagonal,  and  octagonal  shape,  having  the  same 
area  as  the  circular  disk,  were  employed  in  a  limited  number  of  experi¬ 
ments  and  yielded  results  that  differed  imperceptibly  from  those  obtained 
with  circular  disks. 

Vortex  Formation 

When  a  bluff  body  moves  linearly  through  a  fluid  with  a  constant 
velocity,  vortex  formation  and  shedding  occurs  in  the  wake  with  a  period 
defined  by: 

Ts  =  d/Stv,  (20) 

where  St  is  a  Strouhal  number  (nearly  constant  over  a  wide  Reynolds  number 
range).  Bodies  moving  with  different  constant  velocities  shed  vortices 
at  different  periods;  however,  for  steady  flow  conditions  v  =  x/t,  the 
displacement  x  of  the  body  during  each  shedding  cycle  is  the  same,  inde¬ 
pendent  of  the  velocity: 

x  =  d/St.  (21) 

For  a  flat  plate  normal  to  the  flow,  the  Strouhal  number  was  reported 
by  Lehnert26  to  be  0.147  for  Reynolds  numbers  above  4000.  Since  the  vox- 
tex  shedding  from  a  plate  occurs  in  alternate  fashion,  first  from  one  side 
then  from  the  other,  the  distance  of  travel  for  the  shedding  of  a  vortex 
pair  is  6.8  diameters.  Therefore,  there  is  a  3.4  diameter  spacing  between 
shed  vortices. 

If  the  flat  plate  were  oscillated  normal  to  its  surface,  the  ampli¬ 
tude  of  oscillation  y  would  be  just  half  the  total  excursion  in  one  direc¬ 
tion.  To  obtain  a  condition  such  that  no  vortices  were  shed  during  the 
cycle  except  at  the  points  of  direction  reversal,  the  amplitude  of  oscil¬ 
lation  must  be  limited  to  y  =  1.7  d. 

Keulegan  and  Carpenter15  employed  a  stationary  flat  plate  in  oscil¬ 
lating  flow  to  measure  the  hydrodynamic  mass  coefficient  in  a  range  of 
Reynolds  numbers,  based  on  maximum  velocity,  between  5,000  and  15,000. 

The  data  were  correlated  by  the  time  modulus  vt/d  using  the  maximum  veloc¬ 
ity  v  =  2iry/T0,  where  T0  is  the  oscillation  period.  Since  the  time  modu¬ 
lus  can  be  written  in  terms  of  v2/ad  [equation  (14)  and  (16)],  the  data 
of  Keulegan  and  Carpenter  are  shown  in  figure  12. 

The  value  of  Cm,  a  hydrodynamic  mass  coefficient  related  to  the  vol¬ 
ume  of  a  cylinder  of  fluid  with  the  plate  as  its  largest  cross  section, 
increases  with  distance  travelled  (y/d  =  v2/ad)  until  the  vertex  shedding 
cycle  coincides  with  the  oscillation  cycle.  The  hydrodynamic  inertia 
begins  to  decrease  beyond  this  point  as  vortex  shedding  occurs  during  the 
oscillation  cycle.  When  the  amplitude  of  oscillation  is  large  enough 
that  numerous  vortices  are  shed  per  cycle,  a  quasi-steady  condition  is 
reached. 
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FIGURE  12  -  Hydrodynamic  Inertia  Coefficient  Versus  Acceleration 
Modules  for  2  Long  Thin  Plate  in  Oscillatory  Flow 

Stanton  and  Marshall  reported  shedding  of  vortex  rings  from  a  disk 
at  uniform  velocity  (N^e  «  200)  with  a  periodicity  corresponding  to  8 
diameters  of  travel.  Using  Hoemer's27  empirical  function  for  Strouhsl 
number  at  Reynolds  numbers  above  103, 

St  =  0.21/CD3/\  (22) 

the  Strouhal  number  for  the  circular  disk  is  0.192,  and  the  distance  for 
vortex  shedding  is  5.2  diameters.  Vortices  shed  from  the  disk  do  not 
alternate  as  in  the  nlate;  therefore,  one  vortex  would  be  shed  every  5.2 
diameters,  or  when  v2/ad  =  y/d  =2.6,  for  the  oscillating  disk. 

It  was  noted  in  the  oscillation  tests  that  when  the  amplitude  of 
disk  motion  was  in  the  proximity  of  3  diameters  (2y  =  6d),  the  disk  be¬ 
came  unstable  and  tended  to  skew  to  one  side.  This  was  attributed  to 
vortex  shedding  during  the  motion  cycle.  The  effect  was  clearly  detri¬ 
mental  to  the  accuracy  and  credibility  of  the  data  beyond  the  critical 
amplitude,  and  the  experiment  was  limited  to  lower  amplitudes.  It  is 
important  that  utilization  of  the  data  of  the  present  investigation  be 
accompanied  by  the  caveat  against  extrapolation. 

Inversen  and  Balent,  reporting  on  linear  motion  of  disks,  show  scat¬ 
tered  data  above  v2/ad  =  5,  although  the  data  below  v2/ad  =  5  are  well 
correlated.  The  discrepancy  is  attributed  to  side-wall  interference  and 
increasing  inaccuracy  in  this  range,  but  vortex  shedding  would  appear  to 
be  a  highly  probable  source  of  the  scatter. 
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When  the  hydrodynamic  inertia  and  drag  effects  are  coabined  into  a 
single  drag  force,  the  vortex  anomaly,  which  occurs  in  the  former  but  not 
in  the  latter,  has  an  effect  of  far  lesser  Magnitude ,  as  Iversen  and 
Balent's  and  Keulegao  and  Carpenter's  data  illi»strate.  Only  the  cylinder 
data  of  Real eg  an  and  Carpenter  are  indicative  «?f  the  vortex  anomaly  in 
the  total  resistance  coefficient. 

Description  of  an  Oscillation  Cycle 

The  oscillatory  notion  of  a  disk  in  a  fluid  is  not  a  sinple  phenom¬ 
enon,  and  it  is  instructive  to  follow  the  disk  through  the  notion  cycle, 
noting  the  various  fluid-body  interactions. 

Mien  the  disk  is  approaching  the  centerline  of  notion,  the  velocity 
is  nearing  zero.  The  drag  force,  therefore,  is  increasing  to  its  peak 
value  as  inertial  resistance  decreases.  The  moving  disk  has  been  accumu¬ 
lating  and  continues  to  accumulate  a  vortex  ring  that  forms  like  a  dough¬ 
nut  in  the  wake,  enlarging  and  elongating  along  the  direction  of  notion 
as  the  disk  travels  through  the  fluid.  The  energy  associated  with  the 
acceleration  of  the  vortex  and  the  surrounding  fluid,  which  is  not  en¬ 
trained  by  the  disk,  but  nevertheless,  is  set  into  notion  by  it,  are, 
probably,  the  principal  sources  of  hydrodynamic  mass.  Because  the  accel¬ 
eration  at  this  point  is  approaching  zero,  so  too  is  the  inertial  effect 
of  the  hydrodynanic  nass.  At  the  instant  when  the  harmonic  velocity  is 
maximum  and  harnonic  acceleration  is  zero,  the  resistance  force  on  the 
disk  is  entirely  drag,  as  in  the  case  of  the  constant  velocity  linear 
notion,  and  the  hydrodynanic  nass,  although  present,  is  undetectable. 

During  that  portion  of  the  notion  cycle  in  which  the  velocity  in¬ 
creases  from  zero  to  its  maximum,  the  Reynolds  number  describing  the  tui¬ 
tion  increases  from  zero  to  sone  maximum  value  indicative  of  both  laminar 
and  turbulent  conditions  being  experienced.  However,  as  this  motion  cycle 
is  just  one  in  a  long  series  of  oscillations,  the  fluid  has  been  disturbed 
by  previous  cycles,  increasing  the  turbulent  content  of  the  motion. 

As  the  disk  passes  the  centerline  of  motion,  the  velocity  begins  to 
decrease  and  acceleration  effects  become  more  important.  During  this 
quarter  cycle,  drag  decreases  to  zero  and  the  hydrodynamic  inertia  in¬ 
creases  toward  a  maximum.  Reynolds  number  decreases  again,  and  the  com¬ 
ments  made  above  apply.  The  vertex  continues  to  grow.  The  condition  for 
vortex  shedding  has  not  been  attained  as  yet,  since  a  distance  of  travel 
(or  circumstance  of  v2/a)  equivalent  to  six  or  eight  disk  diameters  is 
required  before  sufficient  vortex  mass  and  elongation  for  shedding  has 
been  achieved.  If  the  disk  were  to  exceed  the  vortex  shedding  condition, 
the  process  of  that  shedding  would  cause  an  instability  (presuming  the 
disk  is  not  rigidly  maintained  in  orientation)  to  send  the  disk  askew  and 
cause  considerable  deviation  from  sinusoidal  motion. 
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When  the  disk  attains  the  endpoint  of  the  notion  cycle,  velocity  is 
zero  and  the  disk  is  instantaneously  motionless.  The  vortex  in  the  disk’s 
wake,  however,  is  still  m  motion.  It  may  be  presumed  that  during  the 
previous  quarter-cycle,  this  vortex  has  been  pushing  somewhat  on  the  disk, 
as  the  fluid  entrained  in  the  wake  would  lag  the  disk  in  decreasing  its 
velocity,  but  at  the  point  of  zero  velocity  of  the  disk,  the  vortex  im¬ 
pinges  upon  the  disk,  flowing  over  and  around  it.  As  the  disk  begins  to 
reverse  direction  of  travel,  it  must  overcome  the  inertia  of  the  vortex 
formed  in  its  previous  half-cycle  anu  initiate  the  formation  of  a  new 
vortex  on  the  other  side,  in  its  wake  as  it  travels  again  toward  the  cen¬ 
terline  of  motion,  to  repeat  the  events  of  the  preceding  half-cycle.  It 
is  the  flow-reversal  with  the  momentary  vortex-body  opposition  which  most 
distinguishes  the  oscillatory  case  from  the  linear  one.  The  relative  im¬ 
portance  of  this  end-point-of-motion  effect  to  the  cyclical  hydrodynamic 
inertia  of  the  oscillating  body  is  not  clear,  and  a  detailed  analysis  of 
the  forces  on  the  disk  would  be  required  to  ascertain  this. 

CYLINDERS  AND  TANDEM  DISKS 

The  emphasis  of  the  present  investigation  was  upon  the  determination 
of  the  hydrodynamic  effects  of  harmonically  accelerated  disks;  however, 
cylinders  and  disks  in  tandem  were  also  considered,  employing  the  same 
procedures  heretofore  applied  to  disks. 

Cylinders 


Although  numerous  investigations  of  a  cylinder  in  motion  transverse 
to  its  axis  have  been  made,  the  cylinder  in  axial  motion  has  been  ne¬ 
glected.  Two  and  one-half-inch  diameter  cylinders  with  length-to-diameter 
ratios  (i/d)  of  2,  1,  and  0.025  (disks)  were  oscillated  and  values  of 
hydrodynamic  mass  and  drag  determined. 

The  hydrodynamic  inertia  coefficient  based  upon  the  cylinder  volume, 

Cj  =  mb/eOr/^d2*,  (23) 

was  found  to  vary  with  the  correlation  modulus  v2/ad  and  i/d,  as  shown  in 
figure  13.  The  results  of  Sarpkaya3  with  parallelepipeds  in  small  ampli¬ 
tude  oscillation  and  Patton22  with  parallelepipeds  in  large  amplitude  os¬ 
cillation  are  shown  for  comparison.  These  data  show  the  same  trend  of 
increased  hydrodynamic  mass  with  increased  motion  amplitude  as  the  cylin¬ 
der  data  (y/d  =  v2/ad) .  The  drag  coefficient  was  unchanged  between  i/d  * 

1  and  l/d  =  2,  as  shown  in  figure  14  for  several  values  of  v2/ad.  The 
drag  coefficient  for  cylinders  in  constant  velocity  axial  flow  from 
Hoemer27  is  also  shown.  The  drag  coefficient  for  unsteady  flow  approaches 
the  steady  flow  value  as  v2/ad  increases,  as  in  the  case  of  the  disk. 

A  total  resistance  coefficient  C-p  for  cylinders  was  determined  from 
equation  (7),  using  Cj  and  Cp  from  figures  13  and  14,  and  is  shown  in 
figure  15.  Decrease  of  C-p  with  both  v2/ad  and  i/d  occurred. 


NADC-AE-7120 


Tandem  Disks 


TWo  disks  of  equal  diameters  separated  by  a  distance  s  were  used  to 
determine  the  hydrodynamic  effects  on  oscillated  tandem  disks.  The  disks 
used  were  4  and  6  indies  in  diameter  and  they  were  separated  by  a  1-inch 
diameter  cylindrical  rod. 

Hie  values  of  Cj  obtained  were  based  upon  a  characteristic  volume  of 
d3,  as  in  the  case  of  the  single  disk  {equation  (10)).  The  hydrodynamic 
inertia  coefficient  is  shown  in  figure  16  for  various  values  of  separa¬ 
tion-  to-diameter  ratio  (s/d)  and  vvad.  Wien  s/d  is  very  small,  the  data 
of  the  single  disk  apply.  As  separation  increases  slightly,  the  two  disks 
maintain  a  semblance  of  the  single  disk,  as  the  distance  between  the  disks 
is  too  small,  in  relation  to  the  disk  size,  to  allow  much  flow  in  this 
gap.  At  somewhat  more  separation,  a  decrease  in  C.  occurs  as  destructive 
interference  of  the  vortex  trying  to  develop  behind  the  forward  disk  with 
the  rear  disk  takes  place.  This  interference  is  evident  for  the  range 
0.1  <  s/d  <  0.5,  reaching  a  peak  between  0.25  and  0.35,  with  some  depen¬ 
dence  upon  v2/ad.  Beyond  s/d  *  0.5,  the  tandem  disks  begin  to  exceed  the 
hydrodynamic  mass  of  a  single  disk  until  at  some  large  value  of  s/d  (by 
extrapolation  s/d  >2),  the  tandem  disks  perform  as  two  separate  disks. 

Drag  coefficient  and  total  resistance  coefficient  for  tandem  disks 
arc  shown  in  figures  17  and  18,  respectively.  The  characteristics  of 
these  two  figures  are  much  the  same  as  in  figure  *6.  It  is  noteworthy 
that  Hoeroer27  describes  the  drag  of  tandem  disks  in  steady  flow  by  sum¬ 
mation  of  a  constant  Cp  associated  with  the  forward  disk  and  a  Cp  asso¬ 
ciated  with  the  rear  disk  which  is  negative  in  the  range  0  <  s/d  <2.3 
and  positive  thereafter.  The  drag  coefficient  approaches  that  of  two 
separate  disks  beyond  s/d  *  7.  The  negative  drag  on  the  rear  disk  is 
explained  by  the  suction  effect  on  it  caused  by  the  reduced  dynamic  pres¬ 
sure  within  the  wake  of  the  forward  disk.  As  v2/ad  increases  in  figure 
17,  thus  beginning  to  approach  the  case  of  steady  flow,  the  value  of  s/d 
where  Cq  begins  to  exceed  the  drag  coefficient  of  a  single  disk  increases. 
The  inset  in  figure  18,  Iron  Sfillaarth,  et^  al,28  indicate  a  change  in 
vortex  interaction  between  tandem  disks  at  s7d  =  0.35,  which  approximates 
the  condition  noted  in  the  present  investigation. 
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APPENDIX  A 

DETERMINATION  OF  THE  HYDRODYNAMIC  MASS  AND 
DAMPING  CONSTANT  OF  OSCILLATING  BODIES 


The  motion  of  the  oscillating  body  submerged  in  a  fluid  is  assumed 
to  be  described  by  the  linear  single-degree-of-freedom  equation  for  a 
spring-mass  system, 

my  ♦  cy  ♦  ky  =  F  (A-l) 

The  undamped  natural  frequency  of  the  system  is 

<  =Vk7^  (A-2) 


where  n»p  is  the  physical  mass  and  k  the  spring  constant.  In  the  fluid, 
the  "undamped"  natural  frequency  is 

u„  *-v/k7m  =\/k/(mp  ♦  mh7  (A-3) 

where  m^  is  the  hydrodynamic  mass  and  m  is  the  virtual  mass,  m  *  n»p  +  m^. 
The  damped  natural  frequency  in  the  fluid  is 

und  sVk/m  -  (c/2m)2'  -  (c/2m)2  '  (A-4) 

where  c  is  the  dan.ping  constant  in  the  fluid  (internal  spring  damping  is 
negligible).  The  critical  damping  constant  in  the  fluid  is 

cc  =  2  main  =  2 Vk (mp  ♦  m^)'  (A-5) 

Combining  equation  (A-4)  and  equation  (A-5) 

(c/cc)2  =  1  -  Hd/u*)2  (A-6) 

The  magnification  at  a  frequency  w  is  written 

M  =  y/y0  =  |[l  -  (iu/^)2]2  +  [2(w/u>n)(c/ci:)]2|'1/2  (A-7) 


or 


M  *  y/y0  s  k  [(k-mio2)2  ♦  a)2 c2]"1/2 


(A- 8) 


-  A-l 
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At  u  *  i^, 


Mn  *  cJ2c 


(A-9) 


and  at  Hi  ^ 


M„d  *  {  [l‘-  («nd/^l!)!  *  4(«nd/»„)2(c/cc)^'1/i  (A-10) 


and  sub:  tituting  equation  (A-6) 

1/9 

Mnd  =  |(c/cc)4  *  4(c/cc)2  [1  -  (c/cc)2]| 
Solving  for  (c/cc)2, 

(c/cc)2  =  (4  Wl6  -  H/Mnd2)^ 


(A-ll) 


(A-12) 


and  using  equation  (A-9) 


M„2  *  6/(16  *v/l6  -  12Mnd2) 

(A-13) 

Since  and  Mj,  must  converge  as  damping  appraoches 
creases  as  Mnd  increases,  the  sign  in  equation  (A-13) 
Writing 

zero,  i.e.,  Mn  in¬ 
is  negative  (-) . 

Q  =v/16  -  12/Mnd2 

(A-14) 

(c/cc) 2  *  (4  -  Q)/6 

(A- 15) 

Combining  equation  (A-2)  and  equation  (A-3) 

(wnA>n)2  *  n»/n»p  =  1  ♦  (mh/mp) 

(A- 16) 

Using  equation  (A-4) 

(<%/wnd)2  =  (1  ♦  mh/mp)/[l  -  (c/cc)2] 

(A-17) 

mh  =  "V  f1  -  (c/cc)2]  -  l| 

(A- 18) 

-  A-2  - 
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Substituting  equation  (A-1S), 

%  *  «F  1(2  ♦  Q)/6]  (fu^/und)2  -  1) 
and 

c  *  (■pwn2/wnd)  '  Q)  (2  ♦  Q)/»' 
'diere  Q  is  given  by  equation  (A-14). 


(A-19) 


(A -20) 
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APPENDIX  B 

APPLICATION  Of  THE  DATA  TO  SQNOBUOY  SUSPENSIONS 


The  sonobuoy  suspension  system  nay  be  considered  as  a  spring -mass 
system  governed  by  the  linear  single-degree-of-freedom  equation.  As 
night  be  expected,  however,  application  of  the  results  for  the  hydroay- 
nanic  effects  of  oscillating  disks  to  this  equation  are  not  sinple,  and 
several  alternative  approaches  are  feasible. 

/sumption  A  -  The  physical  parameters  of  the  spring-mass  system 
(k,  m^)V  tre  operational  freouency  («),  the  input  amplitude  (y0),  fluid 
density  (p),  and  disk  diameter  (d)  are  known.  The  amplitude  of  notion 
(y)  is  unknown. 

1.  Alternative  One  -  Iterative  Approach.  The  magnification  equa¬ 
tion  can  be  written 


y  *  Yok  [(k-sea2)2  ♦  «2c2)’!/2 


(B-l) 


a.  Set  m  *  *p  ♦  and  use  figures  6  and  9  where 


cv  =  pCdAv2/2  (B-2) 

or  the  equation 

Cx  *  1.2  (y/d)1/2  (B-3) 

CD  *  2.2  (y/d)-1/2  (B-4) 

or 

b.  Set  m  »  «p  and  use  figure  10  or  the  equation 

Of  =  5,2  (y/d)’1/2  (B-5) 

where 

cv  =  pCTAv2/2  (B-6) 


Estimate  a  value  for  y  and  determine  the  corresponding  values  of  nth  and/or 
c  to  substitute  in  equation  (B-l).  If  the  resulting  y  does  not  equal  the 
original  estimate,  alter  the  estimate  and  repeat  the  procedure  until  the  y 
obtained  from  the  equation  and  the  estimate  are  equal. 


-  B-l 
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2.  Alternative  Two  -  Reduced  variables  approach .  Combining  equa¬ 
tions  (B-2),  (B-3),  and  (B-4)  yields 


1.44c  *  o^w 


(B-7) 


***d  equation  (B-l)  can  be  written  relating  y  to  either  %  or  to  c.  Iter¬ 
ation  can  be  used  to  determine  y,  or  the  equation  (B-3)  or  equations  (B-2) 
and  (B-4)  can  be  used  to  write  a  single  equation  in  y.  This  equation  is 

w%(a2-S)y3  ♦  2oo2(np«2-k)y5/2  *  (k  ay*3)  y2  *  yc2k2  (B-B) 


where 


o  *  1.2  pd5/2  (B-9) 

and 

8  *  o2d5/1.44  (B-10) 

This  equation  nay  be  used  with  an  iteration  approach  soaewhat  siapler  than 
the  previous  alternative,  or  solved  nuaerically. 

3.  Alternative  Three  -  Total  resistance  approach  using  equation 
(B-S)  and  (B-6) 

c2  =  yV  (B— 11) 

where 

Y  *  2^04pd^/2  (B-12) 


Substituting  in  equation  (B-l)  yields 

y3  ♦  f(k  -  mpU2)/Ya)2]2  y2  -  (y0k/Y<*>2)2  =  0  (B-13) 

a  cubic  equation  in  the  form 

y3  ♦  A2y2  -  B2  *  0  (B— 14) 

Again,  iteration  may  be  applied  or  numerical  solution  computed.  Knowing 
A  and  B,  the  equation  may  be  graphed 

F(y)  =  y3  ♦  A2y2  -  32  (B-15) 
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Graphing  F(y)  against  y,  the  function  has  a  aaxiaua  at  y  *  -2A2/3  where 
F(y)  *  [ (4A*/27)  -  B2}  and  a  minimum  at  y  *  0  where  F(y)  *  -B2 .  Therefore, 
F(y)  cros^s  the  y-axis  (Ffy)  *  9)  in  three  places,  with  two  negative  and 
one  positive  value  of  y,  if  A6/B2  >  27/4;  F(y)  is  tangent  to  the  y-axis, 
with  two  equal  negative  values  and  one  positive  val t*e  of  y,  if  A6/82  * 

27/4;  and  F(y)  crosses  the  axis  only  once,  with  two  iaaginary  values  and 
one  positive  real  value  of  y,  if  A6/#2  <  27/4.  In  each  case,  there  exists 
one  real  positive  value  of  y  which  satisfies  equation  (B-14).  The  use  of 
this  graphical  aethod  is  illustrated  in  figure  B-l  where  A2  *  2  and  B2  »  3 
have  been  chosen  arbitrarily.  Intersection  of  F(y)  and  the  y-axis  yields 
the  value  of  y  which  satisfies  the  equation. 

The  equation  (B-14)  can  also  be  used  to  generate  a  graph  as  shown  in 
figure  B-2  where  constant  y  curves  are  plot'  - !  for  values  of  A2  and  . 

For  the  above  example  of  A2  *  2  and  82  *  3,  the  value  of  y  obtained  i*. 
shown  to  agree  with  the  value  obtained  by  graphing  F(y). 

Assumption  B  -  The  physical  parameters  of  the  spring-mass  system 
(k,  *p) ,  cne  operation  frequency  («),  the  input  amplitude  (yc),  fluid 
density  (p),  and  the  desired  amplitude  of  motion  (y)  are  known.  The  disk 
diameter  (d)  is  unknown. 

Several  methods,  similar  to  the  above  solutions  for  y  may  be  used, 
determining  the  value  of  d  from  the  values  of  Of,  and  c.  Equation  (B-13) 
can  be  used  directly  since  y  is  known  to  give 

Y2  *  (y0k/w2)2(i/y3)  -  (k  -  np*2)2(i/y)  (B-16) 

and  using  equation  (B-12) 

d5  «  (yok/pw2)  (1/4. 16y3)  -  f(k  -  BpW2)/pu2)2(l/4.16y)  (B-17) 
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FIGURE  B-l  -  Graphical  Solution  of  the  Total  Resistance  Equation 
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-  Constant  Amplitude  Plot  of  the  Total  Resistance  Equation 
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